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Statistical models on infinite graphs may exhibit inhomogeneous thermodynamic behaviour at 
macroscopic scales. This phenomenon is of geometrical origin and may be properly described in 
terms of spectral partitions into subgraphs with well defined spectral dimensions and spectral weights. 
These subgraphs are shown to be thermodynamically homogeneous and effectively decoupled. 



& , I. INTRODUCTION 



The study of model systems without translation invariance is an interesting and complex subject of modern statisti- 
cal mechanics. A very general description of this situation is in terms of statistical models on graphs, that is on generic 
networks formed by sites, where dynamical variables reside, and links connecting pairwise sites whose variables are 
coupled. This is the direct extension of the typical setup valid for crystalline lattices, which are indeed very special, 
homogeneous graphs. 

On the other hand, graphs are not in general homogeneous and the main question is how these inhomogeneities 
affect physical properties and give rise to relevant changes with respect to lattices. While small scale inhomogeneities 
will affect local properties, one expects that only large scale inhomogeneities are relevant for bulk thermodynamic 
properties. Most likely, the latter properties are those that show universal features which depend only on a few global 
parameters, just as in the case of lattices. The study of such universality requires consideration of infinite graphs 
(with certain natural restrictions given below), where the thermodynamic limit is taken. 

The main relevant geometrical parameter affecting universal properties is the spectral dimension d of an infinite 
graph Q |l]-|^] . It generalizes the Euclidean dimension of lattices to arbitrary real values and is naturally defined from 
the infrared behaviour of the spectral density of the Laplacian operator on Q B. An equivalent definition, the one 
adopted in this work, is in terms of average properties of random walks on Q at large times, that is to say of the 
singularities of the Gaussian model on the same graph |^,|3| . 

On the other hand, the spectral dimension of the whole graph Q, by itself turns out not to be sensitive to macroscopic 
inhomogeneities strong enough to give rise to true thermodynamic inhomogeneities. Indeed it may happen that distinct 
macroscopic parts of an infinite graph exhibit distinct thermodynamic behaviour. We shall show below that such parts 
can be characterized in terms of their own spectral dimension, possibly plus a spectral weight, resulting in an effective 
spectral partition of Q. The crucial point is that these parts form subgraphs which are thermodynamically independent, 
O ■ that is to say completely uncoupled as far as thermodynamic properties are concerned. In other words, inhomogeneous 
| thermodynamic behaviour on the same infinite graph necessarily imply effective decoupling. 

•rH , 

II. INFINITE GRAPHS: BASIC DEFINITIONS, MEASURE AND AVERAGES 

C3 i 

A (unoriented) graph Q || is the ordered couple (G,Gl) formed by a countable set G of vertices (or sites, or 
nodes), that we shall generically indicate with small-case Latin letters, i, j, k, . . . , and a set Gl of unoriented links 
(or bonds) which connect pairwise the sites and are therefore naturally denoted by couples — (J,i)- When the 
set G is finite, Q is a finite graph and we shall denote N the number of vertices of Q. A subgraph Q' of Q is a graph 
such that G'CG and G' L C G L . A subgraph is said to be complete if its has all the available links, that is if, given 
the subset of nodes G' , the subset of links G' L is the largest possible one. 
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A path in Q is a sequence of consecutive links {(i, k)(k, h) . . . (n, m)(m,j)}. A graph is said to be connected, if for 
any two points i,j £ G there is always a path joining them. In the following we will consider only connected graphs. 
The graph topology can be algebraically described by its adjacency matrix A with elements 

The Laplacian matrix L on the graph Q has elements: 

Lij — Z{S{j (2-2) 

where Zi = Y)j Ag , the number of nearest neighbours of i, is called the coordination number (or degree) of site i. 
Here we will consider graphs with z max — sup,; 2j < oo. 

One can also consider a generalization of the adjacency matrix, which corresponds to the ferromagnetic and uni- 
formly bounded coupling Jy, with ^ Aij = 1 and sup Jy < oo, inf > 0. The elements of the generalized 
Laplacian matrix then read: 

Ji&ij J ij (2-3) 

where Ji = Y^j Jij- 

Every connected graph Q is endowed with an intrinsic metric generated by the chemical distance rij which is defined 
as the number of links in the shortest path(s) connecting vertices i and j. 

Let us now consider thermodynamic averages on infinite graphs || . The Van Hove sphere S 0>r C Q of centre o £ G 
and radius r is the complete subgraph of Q containing all i e G whose distance from o is < r and all the links of Q 
joining them. We will call N 0<r the number of vertices contained in S ,r- 

In the thermodynamic limit the average [f]a of a real-valued function / on G is: 

\f\e = lim -r^— fi ( 2 - 4 ) 

™ N °>r i£s., r 

This average does not depend on the choice of the origin o £ G provided / is bounded from below and 

lim J^l = (2.5) 

where |95 , r | is the number of the vertices of the sphere S 0tr connected with the rest of the graph ||. Here we shall 
restrict our attention to graphs with this property. 

The measure \A\ of a subset A C G is the average value [x(-A)]g °f its characteristic function Xi(A) defined by 
Xi(A) = 1 if i e A and X<(^) = if i ^ ^4. The measure of a subset of links G' L C Gl is similarly given by: 

N' r 

\G' L \ = lim (2.6) 

r^oo l\ o r 

where N' L r is the number of links of G' L contained in the sphere S Q _ r . Any two nonzero-measure subsets A and B of 
G are said to be equivalent if their symmetric difference has zero measure, that is \A\ = \B\ = \A n B\. For any given 
nonzero-measure subsets A C G we shall denote {A} its equivalence class. Then A is said to be a representative of 
{A}. With the subgraph 5' defined by the ordered double (G , G' L ), we identify the measure of the subgraph as the 
measure \G\ of its points. 

Given a (nonzero-measure) subset A C G, we define the average on A of any real-valued function f on G as: 

[fU = [x(A)f] G (2.7) 

By definition [J]a is a function only of the equivalence classes, that is [/]a = [f](A\- Moreover, quite evidently 
[f]c = [fU + [/]b whenever C = A U B and \A n S| = 0. 

Given a complete subgraph = (M ,Ml), we denote A4 its complement in Q. This is formed by all points that 
do not belong to M and by all links of Ql which connect them. M is therefore a complete subgraph. We call the pair 
(Ai , M) a partition of order two of Q whenever both M and its complement M are nonzero-measure subsets of G. 

We introduce now the important concept of minimal distance D_(A, B) between any pair A, B of nonzero-measure 
subgraphs of Q such that \A D B\ = 0. It is defined as 

D(A,B) =min(n : \An n B\>0) (2.8) 
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where 



An n B = {i E A : dist(i, B) — n} , dist(i, B) = minr^ (2.9) 

For n — 0, fl n reduces to the usual intersection operator. Notice that, while in general the relation Ap\ n B is not 
symmetric in A, B, the minimal distance is symmetric: D(A,B) = D(B,A). In fact, from the boundedness of Zi, it 
can be shown by induction on n that 

\Bn n A\ > {z max )- n \An n B\ (2.10) 

so that 

\An n B\ > 0=^ \Bn n A\ > (2.11) 

implying our assertion. 

Consider now the minimal distance between the two members of a partition of order two. Suppose Z2(7W, ^Vf ) = 
n > 1; then \M R„_M| > =>■ |M n n _i M| > from the boundedness of z. t . This implies that if D(M,M) is 
finite, then D_(-M, M) = 1. In this case we may say that M and M. are densely interlaced, while in the opposite 
case that they are infinitely separated. From the definition of minimal distance, it follows that if two subgraphs A 
and B of Q are infinitely separated, their common frontier d(A,B) (i.e. the links G Gl with i £ A and j G £>) 
is a zero-measure set. Then the two subgraphs can be disconnected by cutting such a zero-measure set of links. 
This relates the property of infinite separability to the simple separability property defined in ra]. Indeed, the two 
definitions coincide. We shall term separable partition a partition (M. ,M) where M. and M. are infinitely separated. 



III. THE GAUSSIAN MODEL: INFRARED BEHAVIOUR AND THE SPECTRAL DIMENSION 



The Gaussian model on Q is defined |2| by assigning a real-valued random variable (f>i to each node i € G and then 
prescribing the following probability measure 



d^r[4>] = ^- CX P 
Zj t 



i,jeS a ,r 



n ■ 

ies , r 



(3.1) 



for the collection (j) = {fa ; i G S ^ r }. Here Z r is the proper normalization factor, m > is a free parameter and r) 
is the diagonal matrix with elements rjij = rjiSij with the real numbers r\i positive definite and uniformly bounded 
throughout G (that is < 77 min < rj t < ?7 max , Mi &Q). 

The thermodynamic limit is achieved by letting r — > oo and defines a Gaussian measure over the entire 4> — {4>i i * G 
G} which does not depend on the centre of the Van Hove sphere o J(| . The covariance of this Gaussian process reads 



(fafa) = C tj (m 2 ) = (L + m 2 r])^ 
and hence it satisfies by definition the Schwinger-Dyson (SD) equation 

(J, + m 2 r\ i )C li {m 2 ) - 2J JtkC k j(m 2 ) = 8i 



Setting 



13 J 4 + m 2 f]i 



k£G 



Wa = 



J 3 + m 2 ^- 

one obtains the standard connection with the random walk (RW) over Q M : 

oo 



(3.2) 
(3.3) 

(3.4) 
(3.5) 



where the last sum runs over all paths from j to i, each weighted by the product along the path of the one-step 
probabilities in W: 



7 = {i,h- 



,k 2 M,i)=*w[i\ = w ikt _ 1 w kt _ lkt _ 



,w k2kl w klJ 



(3.6) 
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Notice that, as long as m > 0, we have ^2i(W )ij < 1 for any t, namely the walker has a nonzero death probability. 
This implies that CV, is a smooth functions of m 2 for m > e > 0. In the limit m — > the walker never dies and the 



sum over paths in eq. (3.5) is dominated by the infinitely long paths which sample the large scale structure of the 
entire graph ("large scale" refers here to the metric induced by the chemical distance alone). This typically reflects 
itself into a singularity of CV, at m = whose nature does not depend on the detailed form of J,j or rji, as long these 
stay uniformly positive and bounded. 

Of particular importance is the leading singular infrared behaviour, as m 2 — > 0, of the average [C(m 2 )]c of Cu(m 2 ), 
which is a positive definite quantity, over all points i of the graph Q, which we may write in general as 

Smg[C(m 2 )] G ~ cim 2 )*' 2 - 1 (3.7) 

The parameter d is called the spectral dimension of the graph Q and on regular lattices i t co incides with the usual 



Euclidean dimension. Henceforth we shall call spectral weight the coefficient c in eq. (3.7). The name spectral 
dimension is related to the behaviour of the spectral density p(l) of low-lying eigenvalues of the Laplacian L; indeed 
it can be shown that p(l) scales as a power of I for I — * 0, that is p(l) ~ I "v 2-1 . 



IV. LARGE SCALE INHOMOGENEITY: HOMOGENEITY CLASSES AND SPECTRAL CLASSES 

In the study of statistical models one often has to deal with the average [C(m 2 )]A of Cairn r) over a generic positive 
measure subset A C G and in particular one has to consider the leading singular behaviour of [C^m 2 )]^ as m 2 — > 0. 
On regular lattices this singular behaviour is independent of A and it actually coincides with that obtained averaging 
over all points of Q: 

Sing [C{m 2 )] A = Sing [C{m 2 )] G , V4cG, \A\ > (4.1) 

This property arises from the large scale homogeneity of regular lattices due to translation invariance. On graphs, 
where translation invariance is lost, this property can still hold if the inhomogeneity is limited to fi nite scales. More 



generally it may happen that inhomogeneity extends to large scales and the singular parts of eq. (4.1) are different 
on different subsets. However we will prove that such subsets must satisfy very strong topological constraints: a 
large-scale inhomogeneous graph always consists of homogeneous parts joined together by a zero-measure set of links. 
Therefore the splitting of infrared behaviour always corresponds to a macroscopically evident inhomogeneity of the 
graph. 

In this section we will give a rigorous formulation of these statements through the following steps. 

• Let us suppose that the graph Q has indeed a large-scale inhomogeneity that manifests itself through the 
existence of at least one nonzero-measure subset A C G such that, as m 2 — * 0, 

Sing [C(m 2 )] A ~ CAim 2 )^' 2 - 1 (4.2) 

with dA 7^ d. 

• We then define M C G to be a maximally homogeneous (or more briefly maximal) subset with respect to dA 
whenever 

1. \M n A\ > 

2. Sing [C{m 2 )] M ~ c M (m 2 )^'^\ with d M = d A . 

3. For any nonzero-measure subset B C M we have ds = dM- 

4. There exists no B D M such that ds = d,M and \B\ > \M\. 

By this definition it follows that the set of all maximal subsets with respect to dM coincides with the equivalence 
class {M} and we will call it the homogeneity class of dM- 

• Next we prove the 

Theorem 1: The subgraphs M. and its complement M. are infinitely separated, i.e. their minimal distance 
D.{-M,M) is infinite and they define a separable partition of Q. Since this separability is induced by the spectral 
properties embodied by the spectral dimension, we call this a spectral partition (of order two) of Q. 
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• Finally we consider a Gaussian model on the graph M. showing that, from the infinite separability of M and 
M. the spectral dimension of M. is djtf. Therefore, djvf is a property of the graph M. and defines a spectral 
class. This chain of arguments may now be applied to J\A, splitting off a new spectral class if M. has a large 
scale inhomogeneity of the type given above. The process can be repeated until necessary, yielding a complete 
spectral partition of the original graph Q into spectral classes. 

Proof of Theorem 1: 

Let us suppose ad absurdum that D_(M., M) = 1 and therefore that there exists a nonzero-measure subset M' C M 
such that D_(-M , M ') = 1- From the maximality of M it follows that du ^ &m'- Let us consider the random walk 
representation ( |3.5[ ) of Cu(m 2 ) with i 6 M'\ 

Cuim*) = £ W H (4.3) 

7: 2< — i 

Next consider a site k € M whose distance from i is 1. This site exists from the hypothesis D.{M, M') = 1. Then, 
from the sum over paths in the left hand side of (^^) let us retain only the paths containing k. Then, from the 
boundedness and positivity of Jy and rji one gets: 



Cu(m 2 ) > . Ck l ( ™2 (4.4) 
Averaging over M and then over M' we get: 

[C(m 2 )] M , >K[C(m 2 )] M (4.5) 



where if is a positive constant. Now, taking m? —> and using the asymptotic expression for [C(m 2 )] given in (3.7) 
we obtain 

Since this argument applies equally well with M. and interchanged, one gets: 

( m 2)<W2-l > A ,// ^2^-^/2-1 ( 4 ? ) 

which gives du = contradicting the hypothesis. Therefore D.{M,A4) = 00 and A4 and A4 must be infinitely 
separated. 

The infinite separability of M. and M. implies that the two subgraphs can be disconnected by cutting a zero- 
measure set of links. This very peculiar property implies thermodynamic independence, that is the decoupling, in the 
thermodynamic limit, of a model defined on the whole graph Q into two models defined independently on on Ai and 

This applies in particular to the Gaussian m odel , so that the two averages of Cu(m 2 ) on M. and M. are independent 



quantities, each satisfying a relation like eq. (3.7) with two distinct spectral dimensions. Most importantly, to any 
nonzero-measure subset of M there corresponds by construction the same spectral dimension d of A4. We can say 
then that d is a universal property of M. 

V. SPECTRAL WEIGHTS AND SUBCLASSES OF SPECTRAL CLASSES 

In the singular behaviour of [C(m 2 )], inhomogeneity at large scale can appear also in the coefficient of the leading 



infrared part (3.7). However, following the same steps as the previous section, we will show that once again a splitting 
in the value of the coefficient corresponds to a macroscopic inhomogeneity of the graph and that a macroscopically 
homogeneous graph is indeed characterized by universal d and c. Actually in this case the proof is subtler and requires 
some further mathematical steps. 

We first define the spectral subclasses of a given spectral class by looking at the spectral weight ca, proceeding 
along steps similar to those followed above. 
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• Let us suppose that, for a given graph Q belonging to the spectral class characterized by d, there exists at least 
one nonzero-measure subset A C G such that, as m 2 — > 0, 

Smg[C(m 2 )] A ~ ^(m 2 )^ 2 " 1 (5.1) 



with 7^ c, with c given as in eq. (3.7). 

• Then we say that a nonzero-measure subset M C G, which certainly is maximal w.r.t. d, due to its universality, 
is maximal also w.r.t. c A whenever 

1. |M n A| > 

2. Sing [C(m 2 )] M - c M {m 2 )^ 2 -\ with c M = c A . 

3. For any nonzero-measure subset B C M we have cb — cm- 

4. There exists no B D M such that c B = cm and > |M|. 

By this definition it follows that the set of all maximal subsets with respect to cm coincides with the equivalence 
class {M} and we will call it the homogeneity subclass of spectral weight cm- 

• We then prove the 

Theorem 2: The subgraphs A4 and its complement M. are infinitely separated and define a spectral partition 

of a. 

• Following the same steps as the previous section, we then consider a Gaussian model on the graph M. showing 
that, from the infinite separability of M and M, the coefficient of Sing [C(m 2 )]M is cm- Therefore we can say 
that Cm is a universal property of the graph A4 and defines a spectral subclass separated from the rest. 

Proof of Theorem 2: 

To prove this theorem we first need the following lemma: 

Lemma: Within a given spectral subclass, for any subset A of the subclass, the asymptotic form of [C(iti 2 )]a is 
invariant under pre-averaging over any normalized point distribution with nonzero-measure support. In other 
words, if we define 

2 [aC{m 2 )] A 
C(m ) A,a = n ( 5 - 2 ) 

where a.% > on a subset of A with nonzero measure, then again 

Sing [C(m 2 )ka ~ c A {m 2 f/ 2 - 1 (5.3) 

with no dependence at all for c A and d on the distribution a = {at ; i S A}. The proof of this statement is 
elementary: we define the quantities 

/ l = (m 2 )- J / 2+1 a,(m 2 )- CA (5.4) 

Then, by construction, for any e > there exist a S > such that we have \[f] A \ < e as soon as m 2 < 5. Hence we 
also have 

\[af] A \ < (supai)|[/]^| < (supaj)e (5.5) 

ieA ieA 

which immediately implies our assertion. 
Now we can prove Theorem 2: 

Let us suppose ad absurdum that Di-M,M) — 1 and therefore that it exists a nonzero-measure subset M' C M 
such that D_(M. 1 M!) = 1. From the maximality of M it follows that cm ^ c M ,. 

The following proof is given only for d < 4, owing to brevity and physical requirements. Indeed a real structure has 
necessarily a dimension d < 3; moreover, from a purely theoretical point of view, the class of models we have in mind, 
with site variables and link interactions, typically have 4 as an upper critical dimension for the scaling behaviour. 
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Let us consider first the case of a spectral class where [C(m 2 )]a diverges when m 2 — > 0, that is such that d < 2. The 
Schwinger-Dyson equation for Cii[m 2 ] reads: 

(Ji + m\)C u (m 2 ) - J2 J ikC kl {m 2 ) = 1 (5.6) 
keg 



Averaging equation (5J3) over M, we obtain the relation 

[J C] M + m 2 [r, C] m ~ [J ■ C] M = \M\ (5.7) 

where (JC), = JiCu, (rjC)i = T]iCu and (J • C), = J2 k JikCki- We then divide by [J C]m and let m 2 — > 0. Due to 
the divergence of [JC]m we have that, for any e > there exists a (5 > such that, as soon as m < S, 

[J ■ C] M ( rn^ 

Next we set 

Jm',1 = X! Jlk ' (J ' ^)-M"',i = Jik< ^ki (5-9) 
feeM' feeM' 

and use the positivity of Cn — Cik @ to push the above inequality to 

1 e <i [ J m' C 1m , [(J-C)m']m / 5im 

- [JC] M + [jc] m (5 - 10) 



which yields 



Inn K^-%k = 1 (5.11) 

m 2 ^0 [JjQr/ OJm 



Owing to the symmetry of D{Ai, -M'), we may repeat the above steps with M and M' interchanged. Since the 
symmetry of and Cy implies [(J ■ C)m']a/ = [(J ■ C)m]m'j we finally obtain 

m 2 ^0 [Jm OJ M / 

At this stage we apply the lemma given above with a identified with J M ' or Jm , namely 

[J M > G\ M ~ cm [J m >]m {m 2 f 2 - 1 , [J M C} m , ~ c M , [J M ] M ' (m 2 )^ 2 " 1 (5.13) 



But [J M ']m = [Jm\m> so tnat ec l- (5-12) implies cm = c M ,, contradicting our initial hypothesis that D(M.,M.) = 1 
with M maximal. Hence necessarily D_(M,-M) = oo, proving our assertion. 

Let us now consider a spectral class where C(m 2 )a does not diverge in the limit m 2 — > while its first derivative 
with respect to m 2 , C"(to 2 )g, diverges in the same limit. This is the case of a spectral class characterized by a 
spectral dimension 2 < d < 4, where: 



[C'(m 2 )] M , a = [aC r /( ^ )]M ~ -(J/2 - 1) cm (m 2 f 2 - 2 , m 2 -> (5.14) 

,2 



Taking the first derivative with respect to m 2 in the Schwinger-Dyson equation (5.6), we obtain: 

thCu (m 2 ) + m^C^m 2 ) = ]T Jik[C' ki (™ 2 ) - C^(m 2 )] (5.15) 

fees 

which can be averaged over M giving 

[77 C] M + m 2 {v C'] M = [J • C'] M - [J C'] M (5.16) 
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Together with eq. (5.14), this implies 

lim {m 2 ) 2 - d ' 2 {[J ■ C'] M - [J C] M ) = 0+ (5.17) 

mr— >0 

that is, for any e > there exists a S > such that, as soon as m 2 < S 

< £([J • C'] M - [J C'] M ) < e (5.18) 
with £ = (m 2 ) 2 ~ d / 2 . This can be rewritten as 

< [{J ■ C') m ]m ~ [Jm C'] m + [(J • C') M ] M ~ [J M C'] M < C 1 e (5.19) 

Now, since = — rjkCikCkj are the elements of a negative semi-definite matrix, one has that 
[(J • C') m ]m ~ [Jm C'] m > 0. Therefore 

< [(J ■ C') M ] M - [J a , C'] M < r 1 e (5.20) 

Again owing to the symmetry of /^(A^, 7W'), the previous steps can be repeated with M and M interchanged, 
leading to: 

< [(J ■ C') M ] M - [Jm C'] m , < r 1 e (5.21) 
Since [(J • C")j^]m = [(</ • C")m]m, these two relations imply: 

< \[J M > C'] M - [Jm C\ a , | < r 1 e (5-22) 

Eq. ( 5.14 ) entails in the limit m 2 —* 0: 

[Jm> C']m - ~(d/2 - 1) c M [Jm>]m C 1 , [Jm C'] m , ~ -(J/2 - 1) c Al , [J M ]^, T 1 (5.23) 



so that, since [Jm']m = [Jm]m' fr° m ( 5.22 ) one obtains cm — c M >, which contradicts our hypothesis D_(A4, Ai') = 1 
and therefore proves our assertion D(-M, M) = oo. 
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